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Abstract 

We derive the light cone QCD sum rule for the r]NN coupling 
constant g V NN- The contribution from the excited states and the 
continuum is subtracted cleanly through the double Borel transform 
with respect to the two external momenta, p\, p\ = (p — q) 2 . Our 
result is g^nn — (0.3 ± 0.15), which favors small values used in 
literatures. 

PACS Indices: 14.40.Aq; 13.75. Gx; 13.75.Cs 



I. INTRODUCTION 

Nowadays quantum chromodynamics (QCD) is widely believed to be the underlying 
theory of the strong interaction. Yet the non-abelian nature of the gauge group makes 
analytical calculation extremely difficult in the low energy sector. A typical example is 
the various coupling constants of meson nucleon interaction. These couplings are inputs 
for the one boson exchange potentials for the nuclear forces and the analysis of the impor- 
tant pseudoscalar and vector meson photo- and electro-production experiments currently 
underway in MAMI (Mainz) and Spring8 (JHF) etc. For the pion nucleon sector there 
is enough precise data to extract these couplings. Then they are used as inputs to make 
predictions and analyze other experimental data. In the kaon nucleon hyperon sector the 
situation is not so encouraging. But there is still some data available. The worst occurs 
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in the rjNN and rj'NN sector, where knowledge of them is rather poor. In the present 
paper we shall focus on the calculation of rjNN coupling constant. 

There were some theoretical papers on this issue. But the results from various ap- 
proaches differed greatly. With SUf(3) symmetry it was found a V NN = = 3.68 
from the analysis of the nucleon nucleon potential [||]. Similar values was obtained in the 
non-relativistic model 0. From the analysis of forward nucleon nucleon scattering using 
the dispersion relation it was found that a V NN < 1, consistent to be zero ||. In (4| the 
author was able to relate the proton matrix element of flavor singlet current in the large 
N c limit to the pseudoscalar meson nucleon coupling constants, leading to o^nn ~ 1-3. 
Typical values of Oi n NN obtained in fits with one boson exchange potentials range from 3 
to 7 since the eta meson does not contribute significantly to the iViV phase shifts and nu- 
clear binding at normal densities []5J. However this coupling is smaller than 1 and can be 
neglected in the full Bonn potential ||. a V NN extracted from the reaction ix^p — > r\n lies 
between 0.6-1.7 [[?]]. An interesting indirect constraint of a^NN comes from the Ti-rj mix- 
ing amplitude generated by NN loops and neutron proton mass difference using hadronic 
models. In order to let this amplitude agree with results from chiral perturbation theory, 
a V NN is required to be in the range 0.32-0.53 |§. Eta meson photo-production did not 
fix 

o: v nn either. In cx v nn was suggested to around 1.0 or 1.4. Yet a recent analysis of 
more precise eta meson photo-production experiments in Mainz suggested smaller value of 



ctrjNN [10]. In other words, the eta nucleon coupling constant is still very controversial. To 
derive it within an independent and reliable theoretical framework shall prove valuable. 
We shall use the now well developed light cone QCD sum rules (LCQSR) technique to 
calculate a v ^N in this work. Note our approach differs from all the above ones in that it 
starts microscopically from the QCD Lagrangian. 

QCD sum rules (QSR) |]TT[ are successful when applied to the low-lying hadron masses 
and couplings. In this approach the nonperturbative effects are introduced via various 
condensates in the vacuum. The light cone QCD sum rule differs from the conventional 
short-distance QSR in that it is based on the expansion over the twists of the operators. 
The main contribution comes from the lowest twist operator. Matrix elements of nonlocal 
operators sandwiched between a hadronic state and the vacuum defines the hadron wave 
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functions. When the LCQSR is used to calculate the coupling constant, the double Borel 
transformation is always invoked so that the excited states and the continuum contribution 
can be subtracted quite cleanly. Moreover, the final sum rule depends only on the value 
of the hadron wave function at the middle point uq = 1/2 for the diagonal case, which 
is much better known than the whole wave function fl2|l . In the present case our sum 
rules involve with the eta wave function (EWF) <f v {uo = \) etc. These parameters are 
universal in all processes at a given scale. 



We have used QCD sum rules to study the meson nucleon strong interactions. In |U| 
the pion is treated as the external field to analyze the possible isospin symmetry violations 
of the pion nucleon coupling constant. Later the light cone QSR (LCQSR) was employed 
to extract the 7riViV(1535) coupling constant, which was found to be strongly suppressed 
131. The same formalism was extended to the case of vector meson nucleon interaction 



| TEH . The values of the vector and tensor coupling constants and their ratios of pNN 
and ujNN interaction from the LCQSR agree well with the ones from the experimental 
data and the dispersion relation analysis. With the advent of the eta meson distribution 
amplitudes up to twist four [16(1, we are now able to calculate the r/NN coupling constant 
with a theoretically well developed formalism. Although rj meson is a Goldstone boson, 
its mass is not small in the real world and comparable with the typical hadronic scale 
due to the explicit breaking of SUf(3) flavor symmetry. We have included the eta mass 
correction in our calculation. Moreover the eta meson is an isoscalar, which leads to the 
big difference of the LCQSR for the r)NN coupling constant from that for the 7r°NN 
coupling. We arrive at a V NN — ^ 3 ^ L = (0.3 ± 0.15). The numerically small value is due 
to the cancellation between the leading term and mass correction terms. This point can 
be seen clearly in later sections. 

Our paper is organized as follows: Section | is an introduction. We introduce the two 
point function for the r]NN vertex and saturate it with nucleon intermediate states in 
section [n|. The definitions of the eta wave functions (EWF) are also presented. Numerical 
analysis and a short summary is given in the last section. 
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II. THE LCQSR FOR THE r t NN COUPLING 



We start with the two point function 

U( Pl ,p 2 ,q) = i f d^xe^ (0\Tr) p {x)f] P (0)\v(q)) W 



with pi = p, p 2 = p — q and the Ioffe's nucleon interpolating field |17 



Vp(x) = e abc [u a (x)C^u\x)]^d c (x) , (2) 

MV) = e ab ciu b (yh»Cu aT (y)]d c (y)Yl 5 , (3) 

where a, b, c is the color indices and C = 27270 is the charge conjugation matrix. For the 
neutron interpolating field, u <-> d. 

II(pi,p2, q) has the general form 
n(pi,p 2 , q) = F(pf,pl, q 2 )q-f 5 + F x {p\,p\, q 2 )^ b + F 2 (p\ 1 p\ 1 q 2 )p^ 5 + F z {p\,p 2 2l q 2 )o iiU ^q v 

(4) 

The sum rules derived from the chiral even tensor structure yield better results than 
those from the chiral even ones in the QSR analysis of the nucleon mass [|l7j. We shall 
focus on the tensor structure 575 and study the function F{jp\-,p\-,q 2 } as in the QSR 
analysis of the pion nucleon coupling constant. 

The eta nucleon coupling constant g^NN is defined by the rjN interaction Lagrangian: 

£ v nn = g V NNNi^ 5 r]N. . (5) 

At the phenomenological level the eq.(|l|) can be expressed as: 

Uip^q) = a>^(g 2 ) ( ^_ M ^_ M ^ + - • (6) 

where we include only the tensor structure 75^ only. The ellipse denotes the continuum 
and the single pole excited states to nucleon transition contribution. Ajy is the overlapping 
amplitude of the interpolating current i]n(x) with the nucleon state 

(0\Vn(0)\N( P )) = \ N u N (p) (7) 



4 



Neglecting the four particle component of the eta wave function, the expression for 
F(pl,pl, q 2 ) with the tensor structure at the quark level reads, 



i J ci 4 xe^(0|Tr /p (x)r ?p (0)|r / (g)) = 

2i J d A xe^e abc e a ' b ' c, Tr{ lu CiS Tb "\x)C 1 ,iST\x)}^ 

+4i J rf 4 xe^e a6c e a ' b ' c Yr{ 7 ,«S T l b \x)C'7 M (0| M a (x)M a \0)|r ? (g))}7 5 7^^ c '(x)7,75 (8) 

where iS(x) is the full light quark propagator with both perturbative term and contribu- 
tion from vacuum fields |14[] . 

By the operator expansion on the light-cone the matrix element of the nonlocal opera- 
tors between the vacuum and eta state defines the two and three particle eta wave function. 
In order to simplify the notations we use qT^q to denote (uT^u + dT^d — 2^^) / We 
also introduce F v = where is defined as 



<0|g(0)7 M 75?(0)|»7(g)>=i7i J ?M- ( 9 ) 



Up to twist four the Dirac components of this wave function can be written as |L6 

rl . 1 



< 0\q(0h^ 5 q(x)\r](q) >= if^ j du e lwp> [(p v (u) + — m 2 x 2 A(u)} 



+7;fv m v— t du e- iuqx B(u) + 0(x 4 ) , (10) 
2 1 qx Jo 

< 0|g(0)«7 5 g(x)|0 >= fnUr, C du e~ iuc > x VP (u) , (11) 

Jo 



< 0\q(0)a^ 5 q(x)\0 >= -f^q^Xy - q^x^) / du e luqx (p a (u) , (12) 

o Jo 

< 0\q(0)a a p'j 5 g s G l , u (ux)q(x)\ri(q) >= 

if v l^ v V3[(q^qa9u(S - gvQ.a9yfi) - (<? M <?/3fiW - quQpQiia)] J Poj f 3v (ai)e~ igx{ai+va3) , (13) 

< Q\q(Qhnl59sG a p(vx)q(x)\ri(q) >= 

f v m 2 [q,(g ail - - q a (g p , - ^|)] J Pa^^e^^) 

+f v m 2 ^-(q a x f3 - q p x a ) J Va^a^e'^^ (14) 



and 



< 0|g(0)7^G a/3 (ra)g(x)|r/(g) >= 

-if v m 2 v [qp(g a , q a (g„ - )] / Va^aAe'^^ 

-if n ml-^-{q a xp - q p x a ) J Va i ll (a i )e~ i ^+^ . (15) 

The operator G a /3 is the dual of G a p: G a/ 3 = \e a ps p G Sp ; Va.i is defined as Voti = 
dotida>2da 3 5{l — a± — a<i — 0:3). Due to the choice of the gauge x^A^ix) = 0, the path- 
ordered gauge factor Pexp (ig s Jq 1 dux^A^ux)) has been omitted. 

The EWF <f n (u) is of twist two , ip P (u), ip a (u), and {p 3r} are of twist three, while 
A(u), part of B(u) and all the EWFs appearing in eqs. flliD , (|15|) are of twist four. The 



EWFs if(xi, fi) (/i is the renormalization point) describe the distribution in longitudinal 
momenta inside the eta meson, the parameters Xi {J2i x i = 1) representing the fractions 
of the longitudinal momentum carried by the quark, the antiquark and gluon. 



The normalization and definitions of the various constants can be found in |16| . Some 
of them are Jq 1 du tp v (u) = Jq du <4> a {u) = 1, / Uaiip±(ai) = J Pa^n (aj) = etc. 

Since the steps to derive LCQSRs are very similar to those in [|PI|,|I5| , we present final 
sum rule directly. Interested readers may consult the above papers for details. 



\ 2 ^ 

m N X N g vNN e m 1 



F \a^ a (u ) + %>o)]M 2 / (-^-) 



+ Y2^ F n^v r ]3amp 1 [cp 3v \ - —F^ n r] 3 aI 2 [cp 3TI }M 2 f (^) 
+ ^m^(F; + F*){-I 3 [ n ] - I 3 [<p ± ] - I 5 [0 l{ } - J 6 [^||] + h[0 ± ] + / 6 [^]}M 2 /o(^)} , (16) 

n k 

where f n (x) = 1 — e~ x S it is the factor used to subtract the continuum, s is the 

k=0 

continuum threshold. u = M z^ M n , M 2 = ^+m' 1 -> > ^1 are the Borel parameters, and 
ip' a (u ) = ^ftpHi | M=Mo . In order to make comparison with the sum rule for tt° NN coupling 
constant g^NN-, we have labeled the eta meson decay constant F v with the flavor index. 



The functions h[f^\ etc are defined as: 



"0 



<M«o) = - / duB(u) , (17) 
jo 

r u o rl-ua F(ati, Ot 2 , 1 — a± — a 2 ) 



i 1 F =2« Q / aai / cta 2 7- (1 - 2u + ai — a 2 ) , 

J Jo (1 — ai — a 2 ) 



u — ai Jo 1 — u — a 2 

u o r 1 -'^ F(ai, a 2 , 1 - ai - a 2 ) 

aai / aa 2 7- r 

Jo (! — «! — a 2 l 



-2/_ d ai /_ da 2 ^r 2 : " 1 "" 2; , (19) 



J 3 F = 2u / dai / da 2 F a 1; a 2 , 1 - «i - a 2 1 m °— , 20 

Jo Jo (1 — ai — a 2 ) 

/ 4 [F] = 2«„ f° da, da /'""" 2 ' 1 ""'" 2 ' , (21) 
Jo Jo 1 — Qfi — a 2 

I 5 [F] = 2m / da x / da 2 F(ai, a 2 , 1 — a x - a 2 )- , (22) 

Jo Jo 1 — «i — a 2 

h\F] = 2uq / dai I doLzFioLx, 1 — ai — a3,«3) , (23) 
Jo JO 

I 7 F = 2w { / dai / da 2 } , (24) 

JO Uq — Ol\ JO 1 — Mo ~~ a 2 

where F = ip 3r] , <p h <p ± , (p h <p ± . 

III. DISCUSSION 

Since eta meson is an isoscalar, we have F™ = Fff = F v = Replacing the rj index 
with 7i and F^ by f n in (0), we recover the sum rule for g n NN 0]- Note /" = —f£ = f n . 
In other words, the twist four terms involved with three particle pion wave functions 
vanish due to isospin symmetry. The first term in (|i~6D is the leading twist two term. 
The third term is of twist three and of the same sign as the leading term. The second 
term comes from two particle EWF and the remaining terms all come from three particle 
EWFs. Although they are of twist four except the fourth term, their contribution is greatly 
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enhanced by the factor m 2 in contrast with m 2 in the ttNN coupling case. Moreover they 
are of the opposite sign as the leading twist two and three terms, which leads to strong 
cancellation. In other words, large mass and isoscalar structure of eta meson causes g V NN 
to be much smaller than g n NN- 

The sum rule (|T6|) is symmetric and diagonal, which requires the Borel parameters 
Mf = Mf , i.e, Mo = |. The working interval for analyzing the QCD sum rule (^) is 
0.9GeV 2 < M§ < 1.8GeV 2 , a standard choice for analyzing the various QCD sum rules 
associated with the nucleon. In order to diminish the uncertainty due to A^, we shall 
divide by the Ioffe's mass sum rule for the nucleon: 

W#=^)4MVo(-^) + ^-^. (25, 



The various parameters which we adopt are f v = (0.133 ±0.01) GeV [pj|1 , 773 = 0.013, 
a = -4tt 2 < 0|gg|0 >= 0.67 GeV 3 , ^ = 2.13GeV |§ at the scale « = lGeV, s = 
2.25GeV 2 , m N = 0.938GeV, A^ = 0.026GeV 3 JT7J] . 

At uq = I the values of various eta meson wave functions are: ^(uo) = 1.05, A(uq) = 
4.14, <f) B {u ) = 0, (p„(uo) = 1.44, (pf a (uo) = 0, h[(p 3r] ] = 0, h[ n ] = 0.026, I 2 [tp 3v ] = 
-0.9375, I 2 [ n ] = 0, I 2 [<p ± ] = 0, I 3 [ n ] = 0, I 3 [<p ± ] = 0, h[0 l{ ] = -0.313, h[0 l{ ] = -0.032, 
h[0±\ = 0.0396, J fl [^||] = -0.052, I 6 [<p ± ] = 0.044, I 7 [<p ± ] = at u = \ and « = lGeV. 

The dependence on the Borel parameter M 2 of g V NN are shown in FIG 1 with sq = 
2.35,2.25,2.15 GeV 2 . The final sum rule is stable in the working region of the Borel 
parameter M 2 . We obtain: 

g vNN = (1.7 ±0.3). (26) 

In the above numerical analysis we have used relatively large quark condensate value 
< qq >= —(240 ± 10) 3 MeV 3 , which corresponds to a = 0.67GeV 3 . In the literatures 
another value < qq >= —(225 ± 10) 3 MeV 3 and a = 0.55GeV 3 is also used. Since we are 
not able to know very precisely the quark condensate value, we also present the variation 
of g-qNN with M 2 , Sq with a = 0.55GeV 3 in FIG 2. In this case we have, 

^ = (2.1 ±0.3). (27) 

We have included the uncertainty due to the variation of the continuum threshold 
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and the Borel parameter M 2 in (|2q) and (Wnl - in other words, only the errors arising 
from numerical analysis of the sum rule fll6|) are considered. Other sources of uncer- 
tainty include: (1) the truncation of OPE on the light cone at the twist four operators. 
For example the four particle component of EWF is discarded explicitly; (2) the EWFs 
are estimated with QCD sum rule, which also induces some errors; (3) the continuum 
model used in the subtraction of contribution from the higher resonances and continuum 
spectrum; (4) errors in f v etc. 

With all these uncertainties we arrive at 

a vNN = (0.3 ± 0.15) . (28) 

For the 77, 7/ sector instanton effects might be important. It's well known that a large 
part of 7/ mass comes from the U^(l) anomaly. Through 77 — rf mixing instantons also 
affect eta meson mass and decay constant f n . Fortunately we know from phenomeno logical 



analysis that the 77 — 77' mixing angle is about —20 degrees ||18|| . So for eta meson such 
effects may be not so large as in the rj' channel. Direct instantons favor strongly the scalar 
and pseudoscalar channel and might affect the mass sum rules for the mesons in these 
channels. In our QCD sum rule analysis of eta NN coupling constant we have chosen the 
tensor structure (775. Moreover we have used the experimental values for m ni f v as inputs 
instead of invoking the eta meson mass sum rules to extract them. Hence the possible 
correction from instantons is expected to be relatively small. 

In short summary we have calculated the eta nucleon coupling constant with the light 
cone QCD sum rules. The continuum and the excited states contribution is subtracted 
rather cleanly through the double Borel transformation. Our approach differs from all the 
available methods in the extraction of g n NN and starts from the quark gluon level. So it is 
independent and more reliable to some extent. Our result of a V N^ favors the small value. 
Except the nonrelativistic quark model and fits with one boson exchange potentials, other 
approaches tend to yield small values for a V NN- However in such potentials the eta meson 
was treated as some effective degree of freedom to model other multi-meson correlations. 
Hence the eta meson in these potentials can not be related to the real eta meson seen in 
the photo- or electro-production experiments in a simple way. In other words, the cx^nn 
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in these potentials may be not the same quantity as the coupling we have calculated. 
We hope our extraction of a V NN can be used to analyze future eta meson photo- and 
electro-production experiments. 
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Figure Captions 

FIG 1. The sum rule for g V NN as a function of the Borel parameter M 2 with a = 0.67GeV 3 
and the continuum threshold so = 2.35, 2.25, 2.15GeV 2 . 

FIG 2. The same notations as in FIG 1 except a = 0.55GeV . 
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